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^ ! Abstract 

P-H | The validity of a delta-function approximation for the electromagnetic in- 

teraction of relativistic heavy ions is investigated. The production of e + e~ 
| pairs, with electron capture, is used as a test of the approximation. 
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The production of e + e in peripheral collisions of relativistic heavy ions has attracted a 
great amount of theoretical interest due to its non-perturbative character. The calculations 
are hard to perform and it is common to find substantial differences between the cross 
sections calculated within several approaches |]TJ. A good simplification of the problem has 
been found by Baltz and collaborators ||. They have shown that if one makes a gauge 
transformation in the wavefunction of the form ip = exp {— i\ (r, t)}tp', where 
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the interaction induced by the electromagnetic field of an ultrarelativistic particle is gauge 
transformed to (in our units h = c = m e = 1) 

V{p,z,t) = <p{p,z,t){l-va z )-<P{p = 0,z,t){l-a z /v) , (2) 



where (p(p,z,t) is the Lienard-Wiechert potential at a point r = (p, z), generated by a 
relativistic particle with velocity v = v z and impact parameter b, 
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(p, z, t) = 1 Za \{h - pf + 1 2 (z- vtf] ~ V2 . (3) 



In these expressions 7 = (1 — v 2 ) is the Lorentz contraction factor, and a z is the third 
of the Dirac matrices. 

The second part of eq. (fj) acts as a regularization term of the modified potential. It removes 
the divergence at b = 0. This new potential is very useful since the Lorentz contraction yields 
a delta-function in the longitudinal variables when 7 > 1, and b is not too large. Evidently, 
this is a great simplification since delta-function interactions always lead to a considerable 
decrease of integration steps in perturbative as well as in non-perturbative calculations. 
In the formal derivation of the delta-function interaction was obtained by expanding 
the gauge transformed potential (|2]) into multipoles. Further manipulation of the multipole 
expansion and comparison with numerical calculations have shown that (0) can be expressed 
as 

V (p, z,t) = 5(z- t) Za (1 - a z ) In (b ~ p) . (4) 

We will show that this expression can be obtained in a simpler way. The derivation is useful 
to study the validity of the delta-function approximation. In particular we will test the 
approximation in a solvable problem, namely the production of e + e~ pairs in which the 
electron is captured in an orbit around one of the nuclei (bound-free pairs). 

Using the Bethe-integral |J, the potential (|3]) can be written in the form 

I r g-iq.Ugtq.r 

<p(p,z,t) = Za— d 3 q— , (5) 

z7r J q z — v £ q A z 

where u = b + vt and q = (q 4 , q z ). For relativistic particles we can replace (1 — v a) and 
(1 — a/v) by (1 — a) in the interaction (Q). As shown in ref. this amounts to neglect 
a very small (~ 0{1/^ 2 )) piece of the longitudinal part of the interaction. However, it 
is important to keep the other v factors in their respective places, as they give rise to 
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important combinations of 7 factors in the matrix elements. Moreover, the integral in ([5]) 
diverges logarithmically as v — > 1. 

The exact interaction is then given by 

V (p, z, t) = Za(l- a z ) — J d*q e-^ i^-^-J , (6) 

where the denominator of the integrand in (|5|) has been rewritten in terms of 7. The 
interaction given by eq. (JD is a limit of this integral when we set q 2 /j 2 = 0. It is clear from 
the above equation that neglecting this factor yields the delta-function in (f|). However, to 
emphasize the restrictions on q t and q z let us define 

$ (p, z, t) = ^ Z ' f - = - f d 2 q t \ exp (-iq*.b) [exp (iq*.p) - 1] $ 2 (q t , z, t) , (7) 
Za (1 -a z ) re J qf 

where 

= $ i dq z tT~ V ? 2 = ^e-™"*-*! . (8) 
27r 7 qt + qi/T 2 

Now, using limA^oo (A/2) e~ A ' z ' = 5 (x), we see that for 7 — ► 00, $ 2 does not depend on q t , 
and assumes the form of a delta function: <3> 2 (2, t) = 5 (z — vt). 

In this limit, we can write ([?D as 

$(p,z,t)=6(z-t)$ p {p,t) , (9) 



with 



1 /* 1 f dc^£ 

$ p (p, t) = - d 2 q t — exp (-zqt.b) [exp (zq t .p) - 1] = 2 —{Jo [q t \p ~ 6|] - Jo (&&)} 
n J q t J qt 
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where Jo is the cylindrical Bessel function. The integral over each Bessel function diverges, 
but their difference does not. To show this we regularize the integrals by using 
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where Kq is the modified cylindrical Bessel function. Taking the limit k — > 0, and using 
K (ak) ~ In (ak), for small values of ak, we get 



This is the solution of the Coulomb potential of a unit charge in 2-dimensions. An easy way 
to see this is to use Gauss law for the electric field in two dimensions. One obtains E ~ 1/b, 
where b is the distance to the charge. Since E = —d^/db, the logarithmic form of $ is 
evident. 

The above derivation illustrates the validity of the approximation in terms of the transverse 
momentum transfer q t . It should fail for very soft processes, i.e., those for which q t — > 0. 
Also, it requires that q z is small compared to 7. As shown in ref. ||, q z values in the range of 
one up to 7 units of the electron mass contribute appreciably to the integrals involved in the 
production of free, and of bound-free, e + e~ pairs. It is thus important to check the validity 
of the approximation (^) in a concrete case. We will do this for the production of bound-free 
pairs. The full calculation uses the interaction given by equation (H). For comparison, a 
similar calculation with the term £ = q z j^i replaced by zero in the denominator of (|6]) is 
equivalent to the use of the interaction (§) . 

In ref. it was shown that the Coulomb distortion of the positron wavefunction is an 
important effect in calculations of bound-free pair production. The right magnitude of the 
differential cross sections depends on this effect. This is not relevant in our case, since we 
are only interested in the relative change of the cross sections and probabilities by using the 
exact and the delta-function interaction, respectively. Therefore, for simplicity, we will use 
plane-waves for the positron wavefunction and first-order perturbation theory. The energy 
transfer from the field to the created pair is given by u = e + 1, where e is the positron 
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energy, neglecting the atomic binding energy of the captured electron. Using the interaction 
in the form (H), the amplitude for bound-free pair production is given by 



a fi = i dte 1 "* {ip e+ \V\iP e -} 



where the integral over time yields 
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F (Qt) = J d 3 r exp (i^-) [exp (iq f .p) - 1] (ip e + |(1 - a z )\ip e -) 



and £ is now given by £ = lu/jv. 
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Using plane waves for the positron wavefunction and a hydrogenic K-orbital function for the 
electron, the above matrix element is given by 



^(q*) = /(qt)v(l-aju, 



(15) 



where 
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l/<+|Q-p| 



l/a? H + |Q - p| 
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In these equations a# = 1/a = 5.29 x 10 4 fm is the Bohr radius, v (u) is the positron (elec- 
tron) spinor, p = y/ e 2 — 1 is the positron momentum, Q = (q t , u/v), and Qo = (0, uj/v). 

Integrating the square modulus of (TT3] ) over b yields a delta function 5 (q t — qj). Further- 
more, performing the spin averages we get for the differential cross section in terms of the 
positron energy 
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The integral over the positron scattering angle can be done analytically, as well as the 
remaining integral over q t . 

To test the delta-function interaction we define the function 



A (e) = (da/de)t =0 — da/de / (da/de) 
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which does not depend on Z. 

In figure 1 we plot the function A (e) for SPS, RHIC and LHC heavy ion energies. The 
positron energies are given in MeV units. For SPS the above formulation applies directly, 
assuming that the electron is captured in the target and neglecting the atomic screening 
effects. For RHIC and LHC a transformation of ( |T7| ) back to the laboratory system was per- 
formed. We notice that the delta-function interaction works very well for positron energies 
of the order of MeV for SPS and RHIC and up to 100 MeV for LHC. The approximation 
worsens abruptly at a certain positron energy. The value of e where this occurs is a function 
of 7. In fact, we expect that (da/de)^ =0 starts to differ substantially from da/de for £ of 
the order of one, i.e, for e/7 ~ 1. It is thus more appropriate to plot A as a function of 
e/7. This is shown in figure 2 for the same laboratory energies as before. In this figure e/'-f 
is given in units of the electron mass. We see that all curves collapse into approximately a 
single one. The differences between the results for RHIC and for LHC are imperceptible. 
These results show that the calculations with the delta-function interaction differ from the 
calculations with the exact potential for positron energies e > O.I7 mc 2 . For SPS and RHIC 
this implies positron energies of the order of a few MeV, and for LHC a few hundred MeV. 

In the frame of reference of the nucleus where the electron is captured, the positrons move 
in the very forward direction, within an angle of the order of I/7 along the projectile 
incident direction ||J]. Thus, to study the impact parameter dependence of the production 
probabilities we can safely use p = p 2 . The differential probability for pair production is 
given by 
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(19) 



Inserting the term inside brackets of eq. (|16| ) in the integral above and neglecting terms of 
order l/a 2 H , one gets the result 

' * ' f 1 (£&) - V K! ( V b)\ - h -K Q {ib) m , (20) 



where £ = u/'yv and 7] = u/v — p. 

We define another function A (b, e) to test the impact parameter dependence of the delta- 
function interaction, 

A (b, e) = \(dP/dedQ)^ =Q - dP/dedtt] / (dP/dedQ)^ =0 . (21) 

In figure 3 we plot the function A (b,e) for SPS, RHIC and LHC heavy ion energies, as a 
function of 6/7. For comparison we use two values of the positron energy; e = 1 MeV, and 
e = 10 MeV. For e — 1 MeV all curves agree and one observes that the approximation is 
good up to impact parameters of the order of b ~ O.l^y/e. This is confirmed by looking at 
the curves for e = 10 MeV. Then the agreement at SPS bombarding energies is not perfect, 
even for the smaller impact parameters. But, for RHIC and LHC energies the results are 
basically equal. These results originate from the function given in eq. (|20|) which drops 
sharply to zero at £ ~ 1, i.e., at positron energies of the order of b ~ j/e. This is the 
so-called adiabatic limit. The electromagnetic field has photon energy components up to 
where t int is the interaction time. This time equals (6/7) /v ~ 6/7 for relativistic 
collisions |3J- 

In conclusion, the present study has shown that the delta-function interaction yields rea- 
sonable results as long as ub/j < 0.1. As seen in figure 3, this amounts to b < 0.1j/uj. As 
observed in ref. 0, the most effective impact parameters for this process are of the order of 
b ~ 1/m. We also see in figures 1 and 2 that the differential cross sections da/du are well 
described up to energies of the order of O.I7. 

For other situations, e.g., nuclear fragmentation due to the electromagnetic interaction in 
relativistic heavy ion collisions, the most effective impact parameter is given by b ~ R, where 
R ~ 10 fm. We thus expect that the delta-function interaction works well for e i$ O.I7 MeV. 
Note that 7 is the Lorentz factor in the frame of reference of one of the nuclei, i.e., 7 = 2^ — 1, 
where 7 C is the collider Lorentz factor. Thus 7 is huge for RHIC and LHC energies, and the 
approximation works well for all energies of practical interest in nuclear fragmentation. 



The basic idea of the delta-function interaction is that the electromagnetic field of a rel- 
ativistic charge looks like a very thin pancake. Those processes which do not involve too 
large energy transfers, will not be sensitive to the spatial variation of the field. Then the 
delta-function is a good approximation. For typical energy transfers of the order of 10-100 
MeV in nuclear fragmentation, the approximation works well for b < 0.01 — O.I7 fm. To 
calculate total cross sections it is always necessary to account for those large impact param- 
eters at which the delta-function approximation fails. Similar conclusions have been drawn 
in a recent article on projectile-electron loss in reletivistic collisions with atomic targets ||. 
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Figure Captions 

1. Relative difference between the bound-free pair production spectrum calculated with 
the delta-function approximation and with the exact interaction, respectively (see eq. 
(|i~8|)). The comparison is done for SPS, RHIC and LHC heavy ion energies. The 
positron energies are given in MeV units in the laboratory system. 

2. Same as in figure 1, but as a function of the positron energies divided by the Lorentz 
gamma factor. The positron energies are given in units of the electron rest mass. 

3. Relative difference between the bound-free pair production probabilities calculated 
with the delta-function approximation and with the exact interaction, respectively 
(see eq. fl2"T|)). The comparison is done for SPS, RHIC and LHC heavy ion energies 
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and for two different positron energies. The impact parameter divided by the Lorentz 
factor, 6/7, is given in units of the electron Compton wavelength, H/mc. 
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